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ABSTRACT 

Nowadays the number of models aimed at explaining the Type la supernova phenomenon is high 
and discriminating between them is a must-do. In this work we explore the influence of rotation in 
the evolution of the nuclear flame which drives the explosion in the so called gravitational confined 
detonation models. Assuming that the flame starts in a point-likc region slightly above the center of 
the white dwarf (WD) and adding a moderate amount of angular velocity to the star we follow the 
evolution of the deflagration using a smoothed particle hydrodynamics code. We find that the results 
are very dependent on the angle between the rotational axis and the line connecting the initial bubble 
of burned material with the center of the white dwarf at the moment of the ignition. The impact 
of rotation is larger for angles close to 90° because the Coriolis force on a floating element of fluid 
is maximum, and its principal effect is to break the symmetry of the deflagration. Such symmetry 
breaking weakens the convergence of the nuclear flame at the antipodes of the initial ignition volume, 
changing the environmental conditions around the convergence region with respect to non-rotating 
models. These changes seem to disfavor the emergence of a detonation in the compressed volume 
at the antipodes, thus compromising the viability of the so called gravitational confined detonation 
mechanism. 

Subject headings: hydrodynamics - instabilities, rotation - methods -numerical - supernovae: general 


1. INTRODUCTION 

Because of the many connections of Type la super¬ 
novae (SNe la) with fundamental problems in Astro¬ 
physics s uch as, for example, the chemical evolution of 
galaxies (llwamoto et all 199911 or the origi n of the accel¬ 
eration of the universe ( Riess et al.l l!998h . the quest of 
the mechanism behind the explosion is of utmost rele¬ 
vance. There is a general agreement that the supernova 
display involves the explosion of a white dwarf but be¬ 
yond this point there is not an agreement on the details. 
The explosion may be the outcome of the destabilization 
of a white dwarf which approaches the Chandrasekhar- 
mass limit owing to the mass accretion from a nearby 
companion star, a succession of events referred as the 
single degener ated scenario (SD) (IWhelan fc IbenlH97.SI : 
lNomoto| 1198 2b Another route to the explosion involves 
the merging of two white dwarfs with more or less canon¬ 
ical masses Mwn — 0-7 M^ orbiting in a close bi¬ 
nary system (llben fc Tut ukovlll984 iLoren-Aguilar et al.1 
120091: iPakmor et al.ll201.8f l. known as the double degen¬ 
erate scenario (DD), which inclu des the direct c olli sion 
of two WD as a limit ing case dRosswog et al.l l2009al : 
iGarcfa-Senz et all 120131) . Actually, both channels may 
coexist and contribute to th e observed averaged rate of 
SNe la dWang fe Hanl f20l2h . However, each scenario, 
SD or DD, exhibits a considerable degree of degener¬ 
acy because the basic observational properties of SNe 
la are matched by different explosion models. Exclud¬ 
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ing violent mergers and He-detonations, there is not 
a co nsensus about the mechanism driving the explo¬ 
sion ([Branch &; Kh okhlov 1995; Hillcb randt &; Niemeverl 
120001: iHillebrandt et al.l 120131) . Among the different 

mechanisms found in the literature there are either pure 
(subs onic) deflagrations dNomoto et al.lll983: iFink et al.l 
[2013) or deflagrations which at so me critical densit y 
turn into a (supersonic) detonation dKhokhlovl 1 1991 all . 
although the physical agent driving the deflagration to 
detonation transition is uncertain. Another way to det¬ 
onate the WD relies on the gravitational confinement of 
the plasma (the cosmic version of the terrestrial iner¬ 
tial confinement fusion). In these confinement scenarios 
there is a brief deflagrative phase which causes the ex¬ 
pansion of the WD before a large amount of material 
is burnt. The explosion is rather weak which results in 
a failed SNe la. Still the explosion can be revitalized 
in two ways leading to the pulsating and gravitational 
confined models of thermonuclear supernova respectively. 
The pulsating explosion models arise after the fallback 
of some material previously expelled and the formation 
of an accretion shock which compresses and heats the 
core. In steady conditions the ram pressure exerted by 
the infalling material keeps the plasma compressed long 
enough to keep ignited the 12 C + 12 C reaction, which 
under the appropriate physical conditions, may give rise 
to a Chapman-Jouguet detonation which incinerates the 
whole core. Because in these scenarios the explosion 
ensues after a global pulsation of the WD, they were 
called Pulsating Delayed Detonatio n models ( PDD) in 
spherically symm etric calculations (llvanova et al.111974 
iKhokhlovf Il991bf) and Puls ating Reverse Detonations 
(PRD) in three dimensions dBravo fe Garcia- Senzl 120061 : 
iBravo fc Carcfa-Scnz 1 120091 [Bravo et al.ll2009l l. A sec¬ 
ond possibility is that considered in the gravitation al 
confined detonation models (GCD) dPlewa et al.l 12004 
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Plewal 120071: iRopke et all l2007t iTownslev et al.1 120071: 

Jordan et al l 120081 iMeakin et al l [20091 : iSeitenzahl et al l 

20091) whose main features are: 1) the deflagration flame 
ignites in a small volume offset from the center of the 
WD, 2) the bubble of hot fluid accelerates vertically by 
buoyancy in the strong gravitational field of the com¬ 
pact object, 3) soon the bubble made of ashes expands 
over the surface of the star and converges at the opposite 
side from where the initial breakout occurred, 4) as the 
convergence at the antipodes is strong it produces an in¬ 
wardly moving jet which may give rise to the detonation 
in (or close to) the core. 

Until now all GCD models have been calculated as¬ 
suming that the white dwarf is not rotating at the mo¬ 
ment of the explosion. Nevertheless, we know that in 
both, the SD and DD scenarios, the WD could store an 
important amount of rotation because of the angular mo- 
mentum transferred f r om the accretion disc to the WD 
(|Piersanti et al.l 120031 lYoon &; ban veil 12004 I2005D . In 
GCD models the deflagrative phase sets the conditions 
of the detonation at late times, and it is very likely that 
the rotation axis is not perfectly aligned with the sym¬ 
metry plane of the flame. As a consequence, Coriolis 
force will act with different strengths in different regions 
of the burning material. Therefore, to what extent may 
rotation weaken the efficiency of the convergence of the 
ashes at the antipodes? 

Although common proposed scenarios of SNe la im¬ 
ply the rotation of a WD, the number of multidimen¬ 
sional calculations of the explosion which incorporates 
rotation is really scarce . One of them was carried out by 
iSteinmetz et all (|1992D . simulati ng the detonat i on of a 
fast rotating WD. More recently, iPfannes et al.1 (l2010al) 
have calculated the deflagration of a rapidly spinning 
WD in three dimensions, after the central ignition of 
the fuel, and conclude that the yields of the species 
synthesized during the explosion do not match the ob- 
served spectra. In a s ubsequent work the same authors 
(|Pfannes et alJl2~010bD analyzed the detonation of a fast 
rotating WD and conclude that this mechanism could 
explain the existence of some superluminous SNe la. 

In the present work we analyze for the first time the 
possible consequences of including the rotation in the 
GCD scenario. To do that we have carried out several 
3D simulations using a smoothed particle hydrodynam¬ 
ics (SPH) code with moderate resolution. The amount 
of rotation considered in our exploratory study is modest 
U/U c ~ 0.08, where fl c ~ 4.7 s -1 is the critical Keple- 
rian velocity at the equator, to not appreciably modify 
the spherically symmetric density profile of the WD at 
the moment of the explosion. As we will see, even that 
small amount of rotation is able to change the thermody¬ 
namic conditions at the convergence point of the ashes. 
Therefore, future studies dealing with confinement sce¬ 
narios should incorporate rotation as a basic ingredient, 
adding a step of difficulty to an already complex problem. 

In Section 2 we describe the main features of the hy¬ 
drocode, the initial setting and the method to calculate 
the evolution of the nuclear flame (described with more 
detail in the Appendix). We also provide here a brief 
discussion of the profile of the angular velocity adopted 
in the simulations. A detailed description of the hydro- 
dynamic evolution of the four considered models is given 
in Section 3. Finally, Section 4 summarizes the main 


conclusions of our work. 

2. HYDRODYNAMIC METHOD. FLAME HANDLING AND 
INITIAL SETTING 

In the GCD and PRD scenarios, the final explosion of 
the white dwarf takes place after a considerable change 
in the WD radius. Being a Lagrangian method free of 
numerical diffusion, SPH is well suited to ha ndle sys¬ 
tems h ostin g a variet y of dynamical scales (lRosswod 
I2009bl : ISpringell 120101) and to track the complex ge¬ 
ometry of the nuclear flame which powers the explo¬ 
sion. As angular momentum is exactly conserved in 
SPH, this technique is also adequate to handle prob¬ 
lems involving rotation. We are using an updated ver¬ 
sion of the SPH code devised more than one decade ago 
which was applied to simulate the thermonuclear explo¬ 
sion of a W D by different p h ysical mec hanisms either 
deflagrations (IGarcfa-Senz fc Bravol2005D or detonations 
(iBravo fc Garcfa-Senzll2008D . Nevertheless. our new hy¬ 
drocode SPHYNX (Cabezon & Garcfa-Senz, 2015 in 
preparation), incorporates a large number of state-of-the- 
art improvements that are worth to mention. The most 
relevant update concerns gradient estimation, which now 
relies on an integral approach which is more accurate 
than t he traditional metho d in SPH (|Garcfa-Senz et all 
l20H ICabezon et ahl IMl) . To carry out interpola- 
tions SPHYNX mak es use of the sine family of kernels 
(ICabezon et al.ll2008l) . which are more resistant to par¬ 
ticle clustering than the standard cubic spline, therefore 
allowing to increase the number of interpolating particles 
in the SPH summations to reduce the numerical noise. 
The smoothin g leng th h(r,t) is up dated according to the 
method bv iSpringel fc Hernauistl (120021) . which ensures 
both, energy and entropy conservation. In addition the 
code has been parallelized with an hybrid scheme in MPI 
+ openMP, so that we are routinely using ~2x 10 6 par¬ 
ticles. 

The algorithm to track the flame still relies on a 
diffusion-reaction scheme but it incorporates some novel 
improvements. Among them, a phenomenological sub¬ 
grid model to estimate the characteristic turn-over ve¬ 
locity, v t , of the eddies at scales resolved by the hy¬ 
drocode (Appendix A). The effective velocity of the 
flame is Vf = max(vi, Vt) Vi , being the laminar veloc- 
ity of the flame from iTimmes fc Wooslevl (I1992D . At 
the typical conditions considered in our simulations we 
obtain 100 km s -1 < Vf < 500 km s -1 , the former 
being roughly the conductive flame velocity at p ~ 
2 x 10 9 g cm -3 and the latter only achieved in points 
where the shear is high. Different groups approximate 
Vf using a variety of numerical schemes, based either on 
a careful estimatio n of the turb u lent v elocity plus the 
level-set approach dReinecke et al.|[l997i) . or in diffusion- 


reaction schemes either with constant proi 

jagation veloc- 

itv (IGarcfa-Senz & Bravd 20051 Ma et al.l 

2013[) or pow- 

ered by the Rayleigh-Taylor instability 

Gamezo et al.l 


119991 : iPlewa et al.ll26o4 iPlewal[20071) . In any case, the 
gross details of the evolution are insensitive to the pre¬ 
cise value of Vf because the effective flame velocity self- 
regulates through the creation/destruction of the h eat¬ 
excha nging surface between ashes and fuel (iKhokhlovI 
fl995l) . 

The physics included are very simi lar to those re¬ 
cently used by IGarcfa-Senz et all ([2013D to simulate the 
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collision of white dwarfs. The nuclear network is an 
a—chain complemented with carbon and oxygen binary 
reactions. The evolution of the species is calculated im¬ 
plicitly with the temperature to ensure a smooth transi- 
tion to the nuclear-st atistical equilibrium (NSE) regime 
(ICabezon et alj|2004lh Electron captures on protons and 
nuclei have been neglected because their i mpact on the 
dyna mics of the explosion is secondary dMeakin et al.l 
l2009f). Our EOS h as the contributions of electrons 
({Blinnikov et al.lfl99fit l. ions (including Coulomb and po¬ 
larization corrections) and radiation. 

The initial model is a white dwarf with central density 
p c = 2.6xl0 9 g cm -3 , homogeneous temperature T = 2x 
10' K and composition Xq = Xo = 0.5. The integration 
of the Lane-Emden equation for these conditions gives a 
mass Mwd = 1-37 M 0 . The WD was then mapped to a 
three-dimensional configuration of particles and relaxed 
to get rid of the excess of numerical noise before starting 
the hydrodynamic simulations. The maximum resolution 
is ~ 12 km, achieved at the center of the WD. While this 
resolution is high enough to capture the main features 
of the explosion during the deflagration phase it is below 
of what is needed to study the transition to a possible 
detonation at late times. 

All calculations reported in this paper assume that the 
thermonuclear ignition of the WD starts in a single small 
region (a bubble) close to the center of the star and 
propagates to the remaining plasma by hydrodynamic 
instabilities. This initial setting is similar to that cur¬ 
rently used by other groups easing the comp arison of 
results (|Ropke et al.ll2007t iMeakin et al.l 1200911 . On an¬ 
other note, the ignition in a single region at an altitude 
~ 50 km is favored by the re c ent studies bvlZingale et al.l 
(|2011l k lNonaka et alJ (|2012f ): iMalone et ali i 20141 of the 
preignition stage using the MAESTRO code. In several 
of these works slow rotation was shown to affect the long 
conve ctive phas e preceding the ign ition (|Zingale et alJ 
1201 ID . (see also iKuhlen et al.l [200611 . We also assume 
that before the convergence of the ashes the combus¬ 
tion regime is purely deflagrative, that is we have ig¬ 
nored any potential deflagration to detonation transition 
when the flame enters int o the distributed regime at den¬ 
sities p ~ 2 x 10 7 g cm~ :i iKhokhloyfl 991 at [Ciamezo et al.l 

HHI). 


The main goal of this work is to explore the impact 
of rotation on the GCD models. Nevertheless, we have 
restricted the study to moderate rotators to preserve the 
spherical symmetry and equilibrium properties of the ini¬ 
tial model, which is the same for all calculated models. 
Considering faster rotators would need a totally differ¬ 
ent initial setting, being beyond the scope of the present 
manuscript. The main features of the four calculated 
models during the deflagration phase are given in Ta¬ 
ble Q] 

The rotational velocity of a WD at the end of the accre¬ 
tion phase is poorly known. There is some agreement on 
that the core is probably rotating as a rigid body and, be¬ 
yond that, the angular velocity may increase up to a crit¬ 
ical radius R Cl to fi nally decay to Keplerian values at the 
surfac e of the WD (lYoon & Langedl2004 iPfannes et alJ 
l2010aD . Nevertheless in those cases where an efficient 
redistribution of angular mome ntum is assumed th e 
WD may approach rigid rotation (|Piersanti et al.ll2003l h 
Rigid rotators could be the result of having either strong 



Fig. 1.— Dimensionless indicators of the equatorial rotational 
strength for two choices of the angular velocity of the WD. 
Continuum lines are for rigid rotators with Q = 0.4 s~ 1 . Dashed 
lines correspond to f2(r) calculated with Equation [T] 


braking magnetic torques or efficient viscous angular m o- 
mentum transport and/or long time-scales (jPir d l2008h . 

In our exploratory study we have considered four mod¬ 
els (Table [Ij: no rotation, rigid rotators with fL = 
0.4 s^ 1 and Vl x = 0.4 s _1 respectively, and fl x , 


\l x 


n 0 

R c — Rwd 


D 0 = 0.6 

(r — Rwd) 


r < R c , 
r > R c , 


( 1 ) 


where r is the distance to the center and R c = 5 x 10 7 cm. 
This last case is a simple profile consisting of a rigid body 
rotation near the center which decays linearly to a much 
lower value at the surface. The case tt z ^ 0 sets the 
rotating axis aligned to the initial ascending direction of 
the bubble, whereas Q x ^ 0 takes that axis orthogonal 
to the initial displacement of the bubble. It is in the 
latter case where we expect the maximum effect of the 
rotation, mainly driven by the Coriolis force. The value 
of this acceleration, a c or = 2v x f2 acting on a piece 
of a floating bubble is one of the greatest Coriolis forces 
found in Nature. To have an estimate of the order of 
magnitude, we can take D = 0.4 s _1 (models B and C of 
Table[T|) and a typical rising velocity of a bubble v ~ 2.5 x 
10 s cm s _1 . This yields a cor — 2 x 10 s cm s -2 . A huge 
value, almost comparable to the effective acceleration of 
the bubble imparted by the Archimedean force u e // = 



In Figurc[T|we show several magnitudes related to rota¬ 
tion. As we can see, neither the rotational linear velocity 
nor the centripetal force are large enough to significantly 
break the spherical symmetry of the initial model. The 
largest ratio between the linear velocity and the local 
escape velocity is ~ 5%, obtained just on the equator 
at the surface of the WD. The ratio between the cen¬ 
tripetal force and the floating force acting on a rising 
blob g e ff = (A p/p)g(r) ~ 0.2 g(r ), always remains be¬ 
low 3%. Thus, the perturbation introduced by a moder¬ 
ate amount of rotation hardly affects the stability of the 
initial model during the first tenths of second of flame 
propagation. At longer times the nuclear energy released 
during the explosion outpowers the kinetic rotational en¬ 
ergy by several orders of magnitude. The ratio between 
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Fig. 2. — YZ-slice showing a color-map of temperature (in units of 10 9 k) for models A (upper row), B (central row), C (lower row) in 
Table|T]at times t = 0.59 s, /, = 0.97 s and t = 1.18 s respectively. Coordinates are labeled in units of 10 8 cm. 



0 1 2 3 4 5 

t (s) 

Fig. 3. — Evolution of the kinetic energy and the released nuclear 
energy for models A (continuum red line), B (long-dashed line in 
green), C (dashed-blue line), and D (dots in pink). After t = 
2.5 s nuclear reactions were either switched-off or kept switched-on 
(thick lines). 

the rotational kinetic energy and the binding energy of 
the WD at t = 0 is also small, E rot /E^ n < 0.2% in 
all models. However, the ratio between the Coriolis and 
floating forces can be as high as 10% during the first 
tenths of second of flame propagation and should not be 
disregarded. Moreover, the Coriolis force is pushing the 


bubble all the time during its race to the antipodes and 
the cumulative effect may break the symmetry of the 
convergence of the ashes, especially when the rotational 
axis is not aligned with the initial line of displacement of 
the bubble. 

3. HYDRODYNAMIC SIMULATIONS 
3.1. The deflagration phase 

Model A in Table [T| was calculated assuming that the 
WD is not rotating when the first sparks ignite, and 
serves as a reference model. It is also assumed that these 
seminal sparks ignite in a point-like region close to (but 
not at) the center of the core. In all simulated models the 
geometry of the ignited region is almost spherical with 
an initial radius of 20 km. The center of mass of the hot 
bubble is located 60 km above the center of the WD. To 
provoke the ignition the content of the bubble is isochor- 
ically burned to nuclear-statistical equilibrium (NSE), so 
that the temperature jumps to ~ 9 billion degrees, giving 
birth to a thermal wave. We let the flame grow during 
t = 0.25 s keeping the WD structure frozen to ensure 
smooth conditions around the flame. At that time the 
radius of the hot volume has grown up to 40 km. After¬ 
wards, the time is reset to t = 0 s and the hydrodynamic 
calculation of the explosion starts. The bubble expands 
trying to regain equilibrium and its density drops below 
that of the surrounding plasma. It begins to float, ac- 
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Fig. 4.— Same than Figure [21 but at times t = 2.5,3.43,4.25 s respectively. Because the WD has growth in size only the innermost part 
of the object is depicted to highlight the convergence of the ashes towards the antipodes. 


TABLE 1 

Main features of both the initial models and the deflagration phase 


Model 

ZO 

km 

r° b 

km 

^4 x 
s- 1 

ily 

s- 1 

s- 1 

E 

J -'nnc 

10 51 ergs 

M Si 

Mq 

Mjvi 

Mq 

pcimin ) 
10 7 g cm -3 

h c 

km 

PjetiTg — 3) 
10 7 g cm -3 

t(T 9 = 3) 
s 

A 

60 

40 

0.0 

0.0 

0.0 

0.155 

2.5 x 10 -3 

9.0 x 10 -2 

4.6 

49 

1.78 

4.24 

B 

60 

40 

0.0 

0.0 

0.4 

0.157 

2.6 x 10" 3 

9.1 x 10~ 2 

4.5 

49 

1.36 

4.28 

C 

60 

40 

0.4 

0.0 

0.0 

0.155 

3.1 x 10~ 3 

8.9 x 10~ 2 

4.7 

48 

1.16 

4.18 

D 

60 

40 

< 0.6 

0.0 

0.0 

0.139 

2.3 x 10 -3 

8.1 x 10” 2 

5.7 

41 

1.80 

4.01 


Note. — The meaning of of the columns is as follows: zq ignition altitude; r 9 initial radius of the bubble; Q x , Qy, Qz components 
of the angular velocity at t=0 s; E nU c released nuclear energy; M$i , yields of silicon and nickel; p c (mm), h c minimum density 
and resolution achieved at the center of the core of the WD; Pj e t(Tg = 3),t(Tg = 3) density in the colliding region when the 
temperature of a particle made of fuel exceeds 3 x 10 9 K and elapsed time to reach that temperature. 


quiring a buoyant velocity v ~ 2600 km/s at t = 0.5 s. 
Soon the interaction of the rising blob with the environ¬ 
ment moves its spherical geometry to a toroidal one (see 
the leftmost panels in Figure O- From this moment on 
the flame accelerates due to both, the increase of the ef¬ 
fective surface between ashes and fuel, and to the strong 
vorticity induced by the Kclvin-Helmholtz instability at 
the edges of the rising blob. This last effect is handled 
using the toy model described in Appendix A. 

In Figure [3] we show the released nuclear energy, 
E nuc (t) and kinetic energy for all four calculated mod¬ 


els. The evolution of models A, B and C is very similar 
whereas model D produces a bit less nuclear energy. The 
nuclear energy release shows a strong rise at t > 0.8 s and 
begins to saturate at t ~ 1.1 s, owing to the expansion of 
the WD. At the elapsed time t = 1.5 s all nuclear reac¬ 
tions are virtually quenched while the ashes are surfing 
around the core of the WD in their journey to the meet¬ 
ing point at the south pole, FiguresGJQ] At t = 2.5 s the 
nuclear reactions are switched-off from the hydrocode to 
make the analysis of the convergence of the ashes eas¬ 
ier. Although the injected energy E nuc ~ 1.5 x 10 50 ergs 
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T ma x(10 9 K) 


Fig. 5. — Density of the particle as a function of the instantaneous 
maximum value of temperature of fuel particles Tmax achieved in 
the collision region (with nuclear reactions turned off). Model A 
(continuum line in red), model B (dashed line in green), model C 
(short-dashed line in blue) and model D (dots in pink). The elapsed 
time needed to reach Tmax (uppermost curves) can be read on the 
vertical axis at the right. The vertical bars in black shows the 
region where the transition to a detonation may occur. 

is not enough to unbind the star, it provokes a large 
change in its radius, and the central density drops to 
p c = 4.6 x 10' g cm -3 at t = 3.02 s (model A). The value 
of E nuc is in the range of values obtained i n other 3 D 
calculations with simi lar initial conditions (jR.opke et al.1 
EMlFlnk eFIhlEfTl . 

The stages that precede the collision of the ashes for 
models A, B and C are shown in Figure U For the ref¬ 
erence calculation without rotation (model A), the con¬ 
vergence of ashes takes place around the second snap¬ 
shot at time t = 3.43 s. As we can see, despite the 
imprint of hydrodynamical instabilities there is an as¬ 
tonishing symmetry around the line joining the center of 
the remnant with the meeting point. The collision of the 
conical high-velocity stream of gas promotes the birth 
of two jets moving to opposite directions (see the cen¬ 
tral panels in Figure 0 . The basic physics of jets born 
f rom the co llision of fluid streams in both astrophysics 
( Tenorio-Tagle et al.l I1988D and terrestrial laboratories 
( Harlow fe Prachti 119661 ) is well known. Depending on 
the angle of the collision, it may result in one (for low in¬ 
cident angles) or two jets moving in opposite directions. 
Th e latter is jus t what was rep o rted i n the simulations 
by iPlcwal (2007), iMcakin et al.1 (12009H and can also be 
seen in the second and third columns of Figure 0 

According to Table [T] the amount of 56 Ni synthesized 
during the deflagration is around 0.09 M 0 with released 
nuclear energies ~ 0.15 x 10 51 ergs. Thus, in those models 
where the igniting sparks are localized in a tiny, point-like 
region, a transition to a detonation (by the GCD or the 
PRD mechanisms) is necessary to raise the 56 Ni yields 
to standard observed values, ~ 0.6 M 0 . At this point it 
is worth noting that the analysis of the radioactive iso¬ 
topes present in the debris of type la supernova SN2014J 
suggests the existence of a small amo unt of 56 Ni in the 
external layers of the exploding star (iDichl et al.ll2015h . 
Such early ejection could be spread in a truncated coni¬ 
cal ring made of ~ 0.05 M 0 of 56 Vi, detached from the 
main body of the ejecta (Isern et al. 2015, submitted). 
That amount of 56 Ni is curiously close to those shown in 


Table |T} which were produced during the deflagration of 
the WD. Although in our models the ejected 56 Ni is not 
exactly inside a truncated cone its distribution is biased 
towards the hemisphere where the ignition of the WD 
starts. 

The evolution of the rotating models B, C and D pri¬ 
marily depends on the angle between the rotational axis 
and the initial line of flotation of the bubble. When that 
angle is zero the impact in the dynamics is small because 
1) during the first tenths of second the bubble rises ver¬ 
tically, therefore the angle between Vb and f 1 is small 
and so is the Coriolis acceleration, and 2) acoi remains 
axisymmetric during the lateral expansion of the ashes 
preserving the symmetry of the flux. Thus, the evolution 
of model B is similar to model A, as it can be seen in the 
central row of panels in Figures [2] and 01 There are minor 
differences between models A and B, as for example the 
lower value of density achieved in model B at the head of 
the jet when a fiducial temperature, f.e. T = 3 x 10 9 K, 
is crossed (column 12 in Table 0. Also, the elapsed time 
needed to attain for the first time that temperature is 
larger for model B (column 13 in Table [1]). 

Model C has the same initial angular velocity and ki¬ 
netic rotational energy as model B but with fi = T l x go¬ 
ing in an orthogonal direction to the initial displacement 
of the bubble. The overall view of the evolution of model 
C is shown in the third row of Figures [2] and 0 We 
clearly see that the mushroom-like structures are also 
rotating as they are dragged by the substrate. Actually 
they would look similar to that of models A and B in 
the comoving rotation frame. The impact of rotation be¬ 
comes less acute at late times, because the expansion of 
the WD forces f l to drop in order to preserve angular 
momentum, and by t > 1.5 s any angular movement has 
practically ceased while the radial expansion is strong. 
The last two panels of Figure 0 depict the jet forma¬ 
tion for this case. We see that one of the effects of the 
rotation is to blur the component of the jet moving out¬ 
wards, whose direction is now less defined. Interestingly, 
the chunk of the jet moving inwards appears to move ec¬ 
centrically with respect the line defined by the head of 
the jet and the center of the remnant. As commented 
below this behavior is confirmed by a close-up inspection 
around the jet region. Despite of delivering the same 
amount of nuclear energy than in cases A and B model 
C shows some peculiarities at the time of the collision of 
the ashes. According to Figure 0 the density of the fuel 
particle with maximum temperature at the jet region is 
lower than in the non-rotating case. At the fiducial tem¬ 
perature of T max = 3 x 10 9 K there is a 35% reduction 
in density (column 12 in Table [TJ) and the elapsed time 
for achieving T max is a bit shorter. As the conditions to 
form a detonation by the induction mechanism strong ly 
depend on the density of the fuel dNiemever fo Wooslevl 
119971 iDursi fe Timmesll2006lf . such reduction in density 
may have consequences to the further evolution of the 
WD. 

Model D was calculated assuming that il x (r) follows 
the profile given by Equation [l] In this case, the core is 
rotating as a rigid body with Ho = 0.6 s _1 until a ra¬ 
dius of R c = 500 km. Beyond R c a linear decay of Q x 
was assumed until it becomes negligible at the surface of 
the WD. These values of Ho and R c were chosen so that 
the rotational energy of model D at t=0 s is almost the 
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Fig. 6.— Same as Figures l2l and l4l but for model D at times t = 0.6. 1.21 s (upper rows) and t = 2.54,4.21 s (lower rows). 


same as that of models B and C. The evolution of model 
D, summarized in Figure [Gl was not very different than 
that of model C, but in this case the released nuclear 
energy during the deflagration was a slightly lower. As 
a consequence, the value of density in both, the center 
of the WD and the head of the jet, when the temper¬ 
ature Tg = 3 is achieved for first time in the colliding 
region, are larger than in model C (see Table [1}. This 
result supports the common idea that the total amount 
of released nuclear energy is the primary parameter that 
controls the thermodynamical conditions in the colliding 
region, which in turn determines the chances for the det¬ 
onation of the core. Nevertheless, rotation also plays a 
role because, as shown below, the symmetry of the con¬ 
vergence of the ashes is also broken in model D. As in 
model C, the jet is digging tangentially into the core, 
thus changing the environmental conditions at the tip of 
jet. 

Finally, Figure [7] (models A and B) and Figure[8] (mod¬ 
els C and D) show a close-up view of the jet region at 
the time when a detonation could have emerged if nu¬ 
clear reactions were switched-on. In all cases we can 
see the jet traveling towards the core of the remnant. 
In Figure [7] we see that the jet is broad (> 1000 km) 
due to the low resolution in this diluted region. The 
jet is accreting burnt material at its sides through an 
steady accretion shock and the pressure inside the jet is 
roughly balanced by the ram pressure of the infalling 


material. T he detailed st r uctur e of t he jet has been 
descri bed in iMeakin et al.l (|2009D and ISeitenzahl et all 
(12009I) . using a high-resolution axisymmetric approach. 
In models A and B the velocity field inside the jet is 
pointing to the center of the remnant. In these models, 
the jet is digging into the core (see the rightmost pan¬ 
els in Figure O, whereas at the tip of the jet there is 
a compressed region with higher temperature and den¬ 
sity, and is made of fuel. According to previous studies 
it is just in this region where ther e is some chance for a 
detonation flSeitcnzahl et al.ll2009T ) but knowing it unam¬ 
biguously would require a 3D calculation with very high 
resolution (< 500 m) which incorporates many pieces of 
physics (rotation among them). Models C and D rotating 
around a tilted axis show, however, a different behavior. 
According to Figure 0 a jet is also born in the region of 
the collision of the ashes, but now the focusing effect is 
weaker and, more important, the velocity field inside the 
jet is not pointing to the center of the core. This effect 
can also be seen in the color map of density where the jet 
is ’’bitting” the core laterally, which changes the density 
conditions found by precursor shock ahead the tip of the 
jet. Therefore, it seems that even a moderate amount of 
rotation may have an impact on the detonation forma¬ 
tion and future studies should include this ingredient in 
the simulations. 


3.2. Core detonation ? 
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Fig. 7.— Detail of the collision region for model A at time t = 4.229 s (upper panels) and model B at t = 4.289 s (lower panels). Left 
column shows the velocity field superimposed to the temperature (in units of 10® K) color-map. Right column depicts density in units of 
10 7 g cm -3 . 


The initiation of the detonation at the core edge, in¬ 
duced by the inward ly moving jet, h a s bee n analyzed 
with some detail by ISeitenzahl et al.1 (120091 ) assuming 
that the jet is axisymmetric and restricting the calcu¬ 
lation to a small region around the symmetry axis at 
the core interface. Their convergence study, with resolu¬ 
tion between ~ 4 km to ~ 0.1 km, showed that in some 
cases a detonation is formed via the gradient mechanism, 
when a given volume of fuel has a shallow enough den¬ 
sity and temperature distribution so it burns in a time 
lower than the sound crossing time. A detonation occurs 
provided that the size of the burnt zone is larger than 
a critical v alue w hich st rongly depends on density. Ac¬ 
cording to iRonke et al.1 (120071 ) that radius is ~ 120 km 
at p = 3 x 10 6 g cm” 3 and Xg = 2.3, but it decreases 
to ~ 8 km a t p = 10 7 g cm” 3 and T o, = 2.2. In the cal¬ 
culations by ISeitenzahl et al.l (12009D . with a resolution 
of 0.5 km, these criteria are satisfied behind a compres¬ 
sion front that moves transversely to the tip of the jet, 
facilitating the initiation of the detonation. 

Even though our resolution (~ 50 — 70 km, between the 
center of the core and the top of the jet) is not enough 
to elucidate if there is a detonation of the core, we have 
followed the penetration of the jet farther to compare the 
cases with and without rotation. Starting at some point 
at the plateau in Figure [H we resume the calculation 


switching-on the nuclear reactions. In the same figure 
we see that the nuclear energy input becomes again rele¬ 
vant after t= 3.95 s for model A and t= 3.97 s for model 
C. In Figure [9] we show the color-map of temperature 
and density for model A at times t\ = 4.28032 s and 
t 2 = 4.41304 s. It can be seen that the size of the high- 
temperature region grows in time as matter flows into the 
hot volume through the accretion shock. Nevertheless, 
the penetration of the apex of the jet in At = t 2 — t\ is 
hardly perceptible. On the top of the conical jet, but 
detached of it, there is weak shock-wave digging into 
the core. That wave is clearly visible in the density 
maps of Figure [9] especially at t\ = 4.28032 s. A one¬ 
dimensional cut around the symmetry axis showing the 
density and temperature profiles as a function of coordi¬ 
nate Z is shown in Figure [TO] These profiles are moving 
with phase velocity < 1000 km s” 1 , lower than the sound 
speed in the unburnt material, C s ~ 4000 km s” 1 . Thus, 
at our last computed time, t 2 = 4.41304 s, there was not 
a clear signal of detonation in model A. 

In Figure QT] there are depicted the temperature and 
density color-maps of model C at times t\ = 4.28205 s 
and t 2 = 4.40775 s. As in Model A the nuclear reac¬ 
tions were turned-on above t = 2.5 s and the color-map 
focusses on the region around the jet. The head of the 
jet is penetrating obliquely into the core while it grows 
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Fig. 8.— Same as Figurc[7]but for models C and D at times t = 4.277 s and t = 4.21 s. 


in size owing to the nuclear combustion. As in model 
A a precursor shock is formed on top of the jet at time 
t\. This weak shock soon detaches from the subsonic 
jet and moves right into the core, distorting its spherical 
symmetry. Again there was not any clear indication of 
detonation in this region. The evolution of the released 
nuclear energy of models A and C at late times is shown 
in Figure [3] (thick red and blue lines respectively). As 
we can see both curves tends to saturate but the slope 
of the released nuclear energy is lower in model C. This 
is in agreement with the information shown in Figure [5] 
concerning the density of the fuel at a given temperature 
in the collision region. That density is lower in the ro¬ 
tating models and so they are the nuclear reaction rate 
and the released nuclear energy. 

If a detonation does not occur there are two possible 
outcomes for the remnant, 1) It simply remains oscil¬ 
lating while accreting a fraction of the ejected material 
made, for the most part, of unburnt particles. 2) the fall¬ 
back onto the core of some of the previously expulsed 
matter produces and accretion-shock. As the infalling 
material crosses the accretion-shock it is compressed and 
heated, leading to a second chance for a detonation of the 
cor e. This is the PRD ex plosio n mechanism postulated 
by iBravo fe Garcia-Senzl (1200611 . We have tracked the 
evolution of models A and C for a while, until i ~ 12 s. 
In both cases the evolution is quite complex showing sev¬ 
eral episodes of ashes collision while the core is oscillat¬ 


ing. The follow-up of this long phase, probably leading 
to the PRD ignition of the core, is left for a future work. 

4. DISCUSSION AND CONCLUSIONS 

Multidimensional simulations of the thermonuclear ex¬ 
plosion of white dwarfs have, for the most part, assumed 
that the event takes place in a non-rotating progenitor. 
This is actually a wrong hypothesis because the same ac¬ 
cretion process which builds up the mass of the WD also 
transfers angular momentum from the disc to the com¬ 
pact object. As a result the WD could be spinning at an¬ 
gular velociti es as high as fl ~ 2 — 3 s -1 at the moment of 
the ignition dPiersanti et al .1120031 lYoon fe Langedl2005l : 

iDomfnguez et all 2006T) . The reduction in the effective 

gravity caused by the centrifugal force changes the struc¬ 
ture of the progenitor prior the explosion. The den¬ 
sity profile is not longer spherically symmetric and, for 
the same central density, rotating white dwarfs could be 
more massive than their non-rotating counterpart. If the 
amount of rotation is low, to handle the explosion with¬ 
out the intervention of rotation is usually a good ap¬ 
proximation. But in some particular scenarios even a 
slow rotation may have consequences in the development 
of the explosion. This is the case of the GCD scenario 
which relies on the efficiency of the convergence of the 
nuclear flame at the antipodes of the ignition point. A 
good efficiency is achieved provided the deflagration re¬ 
mains axisymmetric with respect to the line defined by 
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Fig. 9.— Same as Figurc[7]but for model A at times t = 4.28032 s and t = 4.41304 s with the nuclear reactions turned-on. 



Fig. 10.— Profiles of temperature and density of model A ob¬ 
tained from a one-dimensional cut around the symmetry line of the 
jet at times t = 4.28032 s and t = 4.41304 s. 

the center of mass of the WD and the ignition point. 
Nonetheless, even a small amount of rotation may break 
the symmetry of the process because the Coriolis force is 
by nature non-axisymmetric and, being proportional to 
the velocity of the hot rising blobs, it can be strong. The 
Coriolis force is also acting in deflagrations arising from 
a multipoint ignition but in that case the stochasticity 
of the ignition may render its final impact less relevant. 


In this work we have studied the impact -via the Corio¬ 
lis force- of having a moderate amount of rotation in the 
fate of the deflagration arising from a point-like ignition. 
Considering angular velocities D ~ 0.4 s -1 is not at odds 
with the hypothesis of spherical symmetry, and allows 
a meaningful comparison of the results between rotating 
and non-rotating models. Such angular velocity is small 
enough to not change the density profile and equilibrium 
properties of the progenitor, but sufficient to appreciably 
push the hot blobs of incinerated material as they rise at 
velocities approaching the local sound speed. A larger 
value of D will make the Coriolis force even stronger 
and, in this sense, our calculations are rather conser¬ 
vative. On another note, moderate spinning WDs could 
be no so rare because: 1) exploding WDs rotating with 
fl < 1 s -1 are expected for compact systems with total 
masses close to the canonical Chandrasekhar-mass limit. 
The physical structure of these WD wou ld be s im ilar to 
that of their non-rotating co unterparts (jPiersanti et al.l 
l2003t iDomineuez et al.l 1200611 . and 2) as commented in 
the introduction the central deflagration of a fast spin¬ 
ning WD leaves a large amount of unburnt matter, whic h 
is incompatible with observations dPfannes et al.l[2010af) . 

The effect of rotation is more pronounced when the 
rotation axis is orthogonal to the ascending line of the 
bubble, reinforcing the idea that the Coriolis force is the 
main agent behind the differences between spinning and 
non spinning models. In this respect, it is worth to men- 
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Same as Figure [7] but for model C at times t = 4.28205 s and t = 4.40775 s with the nuclear reactions turned-on. 


Fig. 11 — 

tion that the value of the angle between the rotation axis 
and the ascending line of the bubble is unknown. In order 
to have a better insight of realistic values of this angle, 
detailed 3D calculations of the progenitor including ro¬ 
tation are needed. 

For a similar released nuclear energy, models with low- 
moderate rotation show two distinctive features: 1) At 
the same elapsed time the temperature and density af¬ 
ter the convergence of the deflagration at the antipodes 
are lower in rotating models, 2) as a result of the col¬ 
lision of the ashes a conical jet is born which grows in 
size through an accretion shock. In non-rotating mod¬ 
els the jet is axisymmetric with its apex heading right 
to the center of the core. In rotating models, however, 
the conical symmetry is not so perfect and the jet is not 
longer pointing to the center of the core. Both features 
above suggest that, along the total released nuclear en¬ 
ergy, rotation is also an important parameter to take into 
account to understand the GCD route to SNe la. 

Several studies have stated that the self-consistent det¬ 
onation of the core of a white dwarf by the GCD mech¬ 
anism is not easy, probably requiring a fine tuning of 
the physical conditions at the convergence region of the 
flame. In particular the released nu clear energy durin g 
the deflagration phase has to be low (iRooke et al.H2007n . 
the geometry of the flame at late times should remain 
as axisym metric as possible to also produce an axisym¬ 
metric jet dMcakin et al.1 120091) and finally, giving birth 


to a detonation requires a precise preconditioning of the 
plasma at the tip of the j et (a shallow gradients of p, T in 
a large enough volume) dScitcnzahl et al][2009 i . Our ex¬ 
ploratory calculations using low-moderate rotators sug¬ 
gest that detonating the core of the WD by the GCD 
mechanism could be even more difficult in spinning mod¬ 
els because the symmetry of the deflagration is broken 
by the Coriolis force. Therefore a deeper understanding 
of the rotational features of a massive WD at the verge 
of the explosion is needed. 
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APPENDIX 

FLAME IMPLEMENTATION. 

Nuclear flames were handled using a reaction-diffusion scheme. Therefore, the actual thickness of the flame was 
artificially enlarged to the spatial local resolution of the simulation, roughly the smoothing-length parameter h. The 
actual conductivity coefficient was conveniently re-scaled to obtain a prescribed velocity for the flame, which was in 
turn computed from a flame subgrid model. Those cases where the scaling factor is one correspond to real microscopic 
laminar flames. The evolution of the flame is basically controlled by the energy equation: 

^7 = 4 ^t+- v '(KVT)+5 ntIC (Al) 

at p z at p 

where k is the thermal conductivity and the remaining symbols have their usual meaning. Now consider the following 
linear transformations: 


r —> r' = a r; t — > t 1 = b t, (A2) 

which maps the coordinate {r, t} framework of the microscopic thermal wave to the laboratory {r' 7 t'} framework of 
the hydrocode. Equation Al can be written as: 

1 1 fjf1 

du=^ P dp + a 2 -V ■ (re V'T)(—) + S nuc {—) (A3) 

p A p 0 0 

On the other hand, for a steady conductive flame moving with velocity v c : 

v “' v= s < A4 > 

Using the scaling transformations A2 in equation A4 gives: 

a v c • V' = b (A5) 

Therefore, only for a = b the scaled flame will move just at the conductive velocity v c . Thus, setting the ratio a/b 
to: 


a = vj_ 
b v c 


(A6) 


leads to a prescribed effective velocity of the flame vp. The practical procedure to implement the flame in the numerical 
scheme is: 


1. Make an estimation of parameter a = (j), where h is the smoothing length an S is the actual microscopic width 
of the flame (jTimmes fc WoosIevlll992( l. 

2. Calculate the effective flame velocity Vf using a subgrid model (see below). The parameter b is then set to 
b = a (^j), where v c is the laminar flame velocity. 

3. Compute the conductivity «;(p, T, Xi) and released nuclear energy rate and evolve the model using equation A3 
as energy equation to describe the evolution of the internal energy and temperature. 

Because the scheme uses both the physical conductivity and nuclear reaction rates calculated through a 14-isotope 
network, it will provide the actual conductive velocity in the limit of very large resolution characterized by a = b ~ 1. 

The nuclear network is able to follow all stages of the combustion, including the nuclear statistical equilibrium 
regime. Nevertheless, in order to save computing time the material is transformed isochorically to NSE elements once 
the temperature has raised over 3.5 billion degrees and density is higher than 5 x 10' g cm -3 . However, once the NSE 
state is achieved its evolution is followed again using the 14-isotope network. Such procedure allows to adequately 
follow the freeze-out of nuclear reactions. 

A subgrid model was built by simply assuming that the size of the velocity fl uctuation at the minimum scale-length 
solvable by the hydrocode is just the effective burning velocity at such scale iDamkohlerl (jl940f ). In order to evaluate 
the characteristic turn-over velocity of eddies with size ~ h we assume an isotropic velocity field around the particle. 
In cylindric coordinates the velocity at a distance r of a given particle is described by: 


v(r, ip) = Kr n r 


(A7) 
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Fig. 12.— Sketch of t he to y-model used to estimate the characteristic turn-over particle velocity. The figure shows the velocity field 
calculated with Equation EH and n = 1 as seen in the comoving frame of a particle located at the coordinate origin. 

where r is the azimuthal unit vector. For n = 1 expression A7 reduces t o rigid rotation (see Figure IT2l) . Cases with 
n = 1/2 and n = 1/3 are representative of RT instability and turbulence (iNicmever fo Wooslcvlll997l ). 

Taking the curl of the above equation in cylindric coordinates and evaluating it at r = h gives: 

Vxc = (n+l)lfr 1 k (A8) 

where k is the unit vector along Z-axis. Constant K in equation A7 becomes: 

K = — h ^"IV x v\ SPH (A9) 

n + 1 

where |V x v\sph stands now for the curl calculated with the hydrocode. Therefore a rough estimation of the turn-over 
velocity at scale h is 


v(h) ~ Kh n = - h IV x v\sph (A10) 

n +1 

The turn-over velocity vt is finally taken Vt = av(h), with 0 < a < 1. In our calculations we use n = 1/3, 
characteristic of turbulence and a = 0.25 fitted so that the total released nuclear energy during the failed explosion 
starting in a non-rotating single bubble is ~ 1.5 x 10 50 ergs. At the moderate angular velocities considered in this 
work the estimated value of v(h) is practically independent of fl. 
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